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Abstract 

The correlation between the pointing direction of ultra high energy cosmic rays and AGN observed by the 
Pierre Auger Collaboration is explained in the framework of acceleration process in AGN. This acceleration process 
is produced by a rotating accretion disk around a black hole that is frozen-in magnetic field. In a result the accretion 
disk is acting as a induction accelerator of cosmic rays. We estimate the resulting magnetic field in the framework of 
the magnetic coupling process. The results of our calculations allow to make the conclusion that the Flat Spectrum 
Radio Quasars appear the effective cosmic accelerators. We estimate also the attenuation of highest-energy cosmic 
rays in a result of their interaction with ambient radiation field. 
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1 Introduction 

Recently, using data collected at the Pierr e Auger observatory during the past 3.7 years, the Pierre Auger Collaboration 
( The Pierre Auger Collaboration . 2007af bl) have demonstrated a correlation between the arrival directions of cosmic 
rays with energy above ~ 6 x lO^^eV and the positions of active galactic nuclei (AGNs) lying at a distance less than 
^ 75 Mpc. A confidence level of anisotropy have been at value of more than 99% through a number of special tests 
( The Pierre Auger Collaboration . 2007b( ). The authors claimed that the observed correlation is compatible with the 
hypothesis that cosmic rays with the highest energies originate from extra-galactic sources located quite close because 
their flux is not significantly attenuated by the Greisen-Zatsepin-Kuz'min effect, i.e. by interaction with the cosmic 
background radiation. Though the present data can't identify AGN as the sources of cosmic rays unambiguously, 
AGNs or objects having a similar spatial distribution can be considered as possible sources. Certainly, these present 
data do not identify as the sources of cosmic rays unambiguousl y. Therefore, other candidate sources which are 
distributed as nearby AGN can't be excluded. Thus. lDermen (j2007l ) has consi dered gamma ray burs ts and only blazar 
AGN as the sources of UHECRs. More original idea has been suggested by iGrib fc Pavlwl (|20n7l ). They suggested 
the hypothesis that dark matter consist of superheavy particles with the mass close to the Grand Unification scale. 
They claimed that some part of these particles, which were created from vacuum by the gravitation of the expanding 
Universe, can be swallowed by the supermassive black holes and can decay on visible particles producing the flow of 
UHECRs observed by the Auger group. 

The interesting look on AGNs as cosmic supercolliders has be en developed by Kardashev and his colleagues 
( Kardashev . 1995t Shatskv . 20031 Ishatskv fc Kardashev , 20021 2006f ). They considered the structure and magnitude 
of electromagnetic field produced by a rotating accretion disk around a black hole. The accretion disk was modeled 
as a torus filled with plasma and the frozen-in magnetic field. It acts as an induction accelerator of cosmic rays along 
the axis of an ac cretion d i sk. 

According to Shatskv ( 2003h the maximum energy of the accelerated charged particle is represented as 



E„ 



(1) 



where R is the distance from a massive black hole to an accretion torus, D, is an angular velocity of a torus, Vg is the 
gravitational radius and fig is the angular velocity at Vg. 

If the basic parameters of the accretion torus near a black hole are of the same order, i.e. f2 ~ f2g and i? ~ r^, the 
Eq.(IT]) transforms to 
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Another p rocess that can be r esponsible for generation of ultra high energy cosmic rays has been most recently 
considered bv lKoide fc Arai They considered the magnetic reconnection as a basic process of energy extraction 

from a rotating supermassive black hole. Magnetic reconnection is one of the important processes for the particle 
acceleration and heating. The relativistic reconnection can be also responsible for generation of ultra high energy 
cosmic rays. Genuinely this process occurs if the magnetic energy exceeds the plasma energy including its rest energy. 
And in this case the magnetic field of AGN takes a key role in the process of acceleration and generation of cosmic 
rays. 

The basic goal of our paper is to obtain the real estimation of the magnetic field strength in the close vicinity of 
accreting supermassive black hole (SMBH) in AGNs and to obtain the more real estimations of Emax from Eq. ([2]) . 



Magnetic coupling process and determination of magnetic field 
strengths of supermassive black holes (SMBH) 



Recently Q (|2002l) : IWang fc Xiao! (|2002[ ): IWang et all (|2003l ): IZhang. Lu fc Zhand (|2005f) . have studied as a possible 
mechanism for transferring energy and angular momentum from rotating black hole to its surrounding disk the magnetic 
coupling process through the closed field lines connecting the Kerr black hol e with its accretion disk. T his process can 
be considered as one of the variants of the Blanford-Znajek (BZ) process (jBlanford fc Znaiekl . Il977l ). It is assumed 
that the disk is stable, perfectly conducting, thin and Keplerian. The magnetic field is assumed to be constant on the 
black hole horizon and to vary as a power law with the radial coordinate of the accretion disk. 

Since the magnetic field on the horizon Bh is brought and held by its surrounding magnetized matter of a disk there 
must exist some relation between the magneti c field strength and accr etion disk and finally the bolometric luminosity 
of AGN. This relation have been obtained by Gan. Wang fc Lil ( 2007t ). It has a form: 
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Here L^oi — sMc^, Rh = 1 + a/I — {a/ MY , M is the accretion rate, e is the radiation conversion efficiency 
(calculated, for example, bv lNovikov fc Thornd (|l973[ ): lKroii^ (|2007l ): [Shapirol (|2007l )). The Eq.® is transformed into 

' Mq \^ /'v\i 1 



Bh = 6.2 X 10'' 
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where {a/MBn) is the Kerr parameter and rj — Lboi/LEdd and Lndd — 1.3 x 1{)'^^{Mbh /Mq) is the Eddington 
luminosity. 

The relation (j4|) can be used for estimation Emax from Eq.([2]). After substituting Eq.(j4]) into Eq.([2|) one gets: 



...... 1.6 X 10^^ 



(5) 



1 - 



3 Estimating maximal energy for specific active galactic nuclei (AGN) 



We are starting with commonly accepted parameters for accreting Kerr supermassive black holes (SMBHs). 

If one chooses for Eq.® the next values of physical parameters: Mbh = IO^Mq, w 1, e w 0.1 and [a/M bh) 
one gets 
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Below we consider more real situation. A syster natic analysis of a large sample of r adio-loud AGN available in the 
BeppoSAX public archive has been performed by iGrandi. Malaguti fc Fiocch I (l2005l ). Their sample includes AGN 



of various types: Narrow Line Radio Galaxies (NLRG), Broad Line Radio Galaxies (BLRG), Steep Spectrum Radio 
Quasars (SSRQ) and Flat Spectrum Radio Quasars (FSRQ) (see Table 6 from their paper). 

Our calculations of the maximal energy magnitude for these types of AGN are presented in Tables 1, 2 and 3. 
One can see from these results that the maximal energy value essentially exceeding the Greisen-Zatsepin-Kuz'min 
(GZK) cutoff is reached only by FSRQ sources. It is well-known that this class is characterized by the presence of 
a non-thermal beamed component. It shows also a featureless continuum and a significantly flatter average spectral 
slope. These facts give evidence that the most effective acceleration process takes place in this type of quasars. 
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Table 1: Magnetic field strength. 



Source 






Bh[G] 


Bh[G] 


Bh[G] 


Bh[G] 


Bh[G] 


Bh[G] 






(t = o.o 

e = 0.057) 


£ = 0.081) 


£ = 0.16) 


(t=0.95 

£ = 0.19) 


= 0.998 

£ = 0.32) 


ifi = 1-0 

£ = 0.42) 


BLRG 


8.64 


-1.59 


lO'' 


9 X 10^ 


8.3 X 10^ 


5.6 X 10^ 


5.3 X 10^ 


5 X 10^ 


NLRG 


8.14 


-2.92 


3.8 X 10^ 


3.4 X 10^ 


3.2 X 10^ 


3.2 X 10^ 


3 X 10^ 


2.8 X 10^ 


SSRQ 


9.29 


-0.90 


10" 


9.4 X 10-'' 


8.7 X 10^ 


8.7 X 10^ 


8.3 X 10^ 


7.7 X 10^ 


FSRQ 


9.01 


0.44 


6.7 X lO'^ 


6.1 X lO'^ 


5.6 X lO'^ 


5.6 X lO'^ 


5.4 X 10"* 


5 X 10"* 


Seifert 1 


7.23 


-0.59 


1.6 X 10^ 


1.4 X 10^ 


1.3 X 10^ 


1.3 X 10^ 


1.3 X 10^ 


1.2 X 10^ 



Table 2: Maximal energy of cosmic rays. 



Source 






Emax[eV] 

(t-0.0 
£ = 0.057) 


Emax[eV] 

(f -0.5 
£ = 0.081) 


Emax[eV] 

= 9 
£ = 0.16) 


Emax[eV] 

(ft = 0.95 
£ = 0.19) 


Emax[eV] 

(ft = 0.998 
£ = 0.32) 


Emax[eV] 

£ = 0.42) 


BLRG 


8.64 


-1.59 


lOi*^ 


9.8 X 10^** 


9 X lO^** 


9 X 10^** 


8.6 X 10^** 


8 X 10^** 


NLRG 


8.14 


-2.92 


1.3 X IQi** 


1.2 X IQi* 


1018 


IQis 


IQis 


9 X 10^^ 


SSRQ 


9.29 


-0.90 


5 X lO^^ 


4.6 X lO^^ 


4.2 X 1019 


4.2 X lO^^ 


4 X lO^^ 


3.8 X lO^^ 


FSRQ 


9.01 


0.44 


1.7 X 10^° 


1.6 X 10^° 


1.4 X 10^° 


1.4 X 10^° 


1.37 X 10^° 


1.3 X 10^° 


Seifert 1 


7.23 


-0.59 


7 X 10^^ 


6 X 10^^ 


5.7 X 10^^ 


5.7 X 10^^ 


5.4 X 10^^ 


5 X 10^^ 



Table 3: Retrograde motion = -0.9, s = 0.039). 



Source 


Bh{G] 


ErnaxleV] 


BLRG 


1.7 X lO'' 


1.8 X lO^y 


NLRG 


6.5 X 10^ 


2.24 X 10^^ 


SSRQ 


1.76 X 10'* 


8.6 X lO^'' 


FSRQ 


1.6 X 10 "' 


4,2 X 10^" 


Sc-ifcrt 1 


2.7 X 10'' 


1,2 X 10^9 
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4 The interaction of ultra high energy cosmic rays with intrinsic radia- 
tion of active galactic nuclei and quasars 

The basic problem for the mechanism considered here is the attenuation of the cosmic ray flux above 5 x lO-'^^eV 
in a resuh of interaction of this flux with the ambient radiation field around of a supermassive black hole. It is well 
known that, namely, this attenuation by extragalactic background radiation produces the so-called GZK effect. The 
similar effect by ambient radiation of AGN or QSO should be estimated in our situation. 

We estimate of the optical thickness r respect to the interaction of the cosmic ray flux with the ambient radiation 
of AGN and QSO. The value of the optical thickness respect to the P7 process is 

(7) 



ATTcRksT 

where L is the radiation luminosity, R is the characteristic scale distance of the radiation field, T is the temperature of 
a radiation field, a i s the cross-section of the interac tion o f cosm ic protons with the radiation field that is ct w lO^^^cm^ 



( Atovan &: Permed (j2003,), see also .Revnoso fc Romerol ( 20081 )) 



The temperature T is determined by the temperature of an accretion disk T,, and it is determined by (jBonning et al 



20071) : 

1/4 / T \ 1/4 / n \ 3/4 
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where i?i„ is the radius of an accretion disk. 

The relation ([5]) corresponds to the standard accretion model of Shakura and Sunyaev. The value of Tq is determined 
at the radius Rin that is closed to the event horizon and corresponds to X-ray temperature, i.e. Tq w lO^K. 

We consider two situations: (a) interaction of the cosmic ray flux with X-ray radiation fleld and (b) with optical 
and UV radiation fleld. 

The radiation flux is falling down with a distance as R^^ . We determine the characteristic scale distance for X-ray 
radiation field as R~ 3Rg, where Rg is the Schwarzschild radius. For the optical and UV-radiation the characteristic 
scale-distance is ~ (10^ 10^)i?g. In a result we obtain: 



12.8( — ) (^) —;R^rRg (9) 



This Eq.® means that r ^ {Mq/MbhY^^ , i.e. the increase of the mass of SMBH decreases the optical thickness. 

For A4bh — 10^A/q and typical values L — = 10'''^erg/s and Lopt.uv ~ 10^'*erg/s we obtain from Eq.®: 
T = Tx ~ 1.2 and r — Topt.uv ~ 6. These values demonstrate that the escape of a part of accelerated protons of UHE 
from the nearest environment of SMBH is very likely. It shoul d be noted that the co smic ray luminosity of the active 



galactic nuclei constitutes a fraction ~ 10 ''of the optical one (jPedenko et al.l . |2008| ) . This result does not contradict 
to our estimations. 



5 Conclusions 



The correlation between the pointing directions of ultra high energy cosmic rays observed by the High Resolution 
Fly's Eye experiment and Active Galactic Nuclei visible from its nort hern hemisphere l o cation can be explained in the 
fram ework of acceleration process in AGN as cosmic supercollider. Kardashev (1995); Shatskv fc Kardashev ( 20021 
l2006l) suggested that this acceleration process produced by a rotating accretion disk around a black hole that is 
frozen-in magnetic fleld. This accretion disk is acting as an induction accelerator of cosmic rays. 

We estimated the resulting magnetic field in the framework of the magnetic coupling process. The results of our 
calculation allow to make the conclusion that only the Flat Spectrum Radio Quasars can be revealed as the effective 
cosmic accelerators. Only for these objects the maximal energy of cosmic rays can essentially exceed the GZK cutoff. 

Our conclusions are not totally exclude another probable scenario of UHECR generation. We mean the process 
of generation of cosmic ray particles with ultra-high energy in a result of accelerat ion in jets and extended envelopes 



surrounding AGN. This scenario was recently developed by a number of authors (jPermer et al.l . l2008t lAllard et al 



20081: iBerezhkol I2OO8 I: iBierman et al.l . l2008l ) . The decisive choice between various probable scenarios will be done in 
the result of future cosmic observations and move details confirmation of the Pierre Auger Collaboration data. 
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